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Abstract 

The multipole response of nuclei at temperatures T = — 2 MeV is studied using 
a self-consistent finite-temperature RPA (random phase approximation) based on 
relativistic energy density functionals. Illustrative calculations are performed for the 
isoscalar monopole and isovector dipole modes and, in particular, the evolution of 
low-energy excitations with temperature is analyzed, including the modification of 
pygmy structures. Both for the monopole and dipole modes, in the temperature 
range T = 1 — 2 MeV additional transition strength appears at low energies because 
of thermal unblocking of single-particle orbitals close to the Fermi level. A concen- 
tration of dipole strength around 10 MeV excitation energy is predicted in 60 ' 62 Ni, 
where no low-energy excitations occur at zero temperature. The principal effect of 
finite temperature on low-energy strength that is already present at zero tempera- 
ture, e.g. in 68 Ni and 132 Sn, is the spreading of this structure to even lower energy 
and the appearance of states that correspond to thermally unblocked transitions. 
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The multipole response of nuclei at energies close to the neutron emission thresh- 
old, and the possible occurrence of exotic low-energy modes of excitation, present a 
rapidly growing field of research [lj]. For instance, in nuclei with a pronounced neu- 
tron excess, the phenomenon of pygmy dipole resonance (PDR), a low-energy mode 
that corresponds to oscillations of a weakly-bound neutron skin against the isospin 
saturated proton-neutron core, has been analyzed in a number of experimental and 
theoretical studies. In particular, recent (7, 7') and (a, a'j) experiments 0, 0, 0] 
have disclosed interesting features of the concentration of low-energy electric-dipole 
excitations. Exotic nuclear dynamics has also been explored in experiments with ra- 
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dioactive ion beams, extending the studies of low-lying dipole response toward nuclei 
with more pronounced neutron excess 0, 0] . 

Features of PDR in medium-heavy and heavy nuclei place constraints on the 
nuclear matter symmetry energy, and can also provide information on the size of 
the neutron skin 0, S). A concentration of electric-dipole strength at low energy is 
interesting not only as a unique structure phenomenon, but it could also play an im- 
portant role in astrophysical processes. For instance, it could affect neutron capture 
rates in r-process nucleosynthesis 0, 10]. Reactions relevant for the r-process take 
place in stellar environment at finite temperature. Calculations of neutron capture 
rates on neutron-rich nuclei have typically included temperature only as a parame- 
ter of the Lorentzian function used to fold the El strength distribution 0, Eil| . It 
would, of course, be useful to perform these calculations using a microscopic frame- 
work that provides a consistent description of nuclear ground-state properties and 
excitations at finite temperature. In particular, temperature induced modifications 
of low-energy dipole strength close to the neutron separation energy, e.g the thermal 
unblocking of additional transitions, should be described in a consistent theoretical 
framework. 

In this work we introduce a fully self-consistent finite-temperature RPA (random 
phase approximation) framework based on relativistic energy density functionals, and 
explore the evolution of multipole response as a function of temperature. The goal 
is to analyze how low-energy excitations evolve with temperature, the modification 
of the structure of PDR, and the possible occurrence of excitation modes that are 
not present in nuclei at zero temperature. 

The multipole response of hot nuclei has been studied using a variety of theoret- 



ical approaches, including e.g. the RPA based on schematic interactions [12], [13|, [14 



15, 16 



Hartree-Fock plus RPA with Skyrme effective inter- 

The 



linear response theory 

actions [17|, RPA models extended with the inclusion of collision terms 
decay of hot nuclei and thermal damping of giant resonances have also been explored 



18 



in numerous studies, e.g. in Refs. [19|, |20|, |2l|, |22], |23|, |24|, |25j, as well as sum rules at 



finite temperature 26j, |27|. Finite temperature Hartree-Fock Bogoliubov (HFB) the- 
ory [i^] has been used as the basis for the corresponding quasiparticle RPA (23, 3^ . 
Experimental studies of giant resonances in hot nuclei are carried out by measuring 
high e nergy 7 -rays from the decay of resonant states built on excited states of target 



nuclei 31, 32, 33, 34 



The self-consistent finite temperature relativistic random phase approximation 
(FTRRPA) is formulated in the single-nucleon basis of the relativistic mean-field 
(RMF) model at finite temperature (FTRMF). This model has been introduced in 



Ref. [35|, based on the nonlinear Lagrangian NL3 [36j. In the theoretical frame- 



2 



work that we use in the present analysis relativistic effective interactions are imple- 
mented self-consistently, i.e., both the FTRMF equations and the matrix equations of 
FTRRPA are based on the same relativistic energy density functional with medium- 
dependent meson-nucleon couplings 37] - Of course, the description of open-shell 



nuclei necessitates a consistent treatment of pairing_correlations like, for instance, in 



the finite temperature HFB+QRPA framework [29f| . However, in nuclei the phase 
transition from a superfluid to a normal state occurs at temperatures T m 0.5 — 1 



MeV |28l. l30t] . whereas for temperatures above T ~ 4 MeV contributions from states 



in the continuum become large, and additional subtraction schemes have to be im 



plemented to remove the contributions of the external nucleon gas 38]. Therefore, 
we consider the temperatures in the range T = 1 — 2 MeV, for which the FTRMF 
plus FTRRPA should provide an accurate description of the evolution of multipole 
response. 

The FTRRPA represents the small amplitude limit of the time-dependent rela- 
tivistic mean-field theory at finite temperature. Starting from the response of the 
time-dependent density matrix pit) to an external field fit), the equation of motion 
for the density operator reads 



2d t p=[h[p} + fit),p}. (1) 

In the small amplitude limit the density matrix is expanded to linear order 

p(t)=p° + Sp(t), (2) 
where p° denotes the stationary ground-state density. 



pli = km 



[1 + expi^-j^)} 1 for states in the Fermi sea (index k, I) 







for unoccupied states in the Dirac sea (index a) 

(3) 

and includes the thermal occupation factors of single-particle states n^. The resulting 
FTRRPA equations read 
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The RRPA matrices A and B are composed of matrix elements of the residual 
interaction V, derived from an effective Lagrangian density with medium-dependent 
meson-nucleon couplings [37|], as well as certain combinations of thermal occupation 
factors rife. The diagonal matrix elements contain differences of single particle ener- 
gies e m — Cj, where e, < e m . Because of finite temperature, the configuration space 
includes particle-hole (ph), particle-particle (pp), and hole-hole (hh) pairs. In addi- 
tion to configurations of positive energy, the FTRRPA configuration space must also 
contain pair-configurations formed from the fully or partially occupied states of posi- 
tive energy and the empty negative-energy states from the Dirac sea. The inclusion of 
configurations built from occupied positive-energy states and empty negative-energy 
states is essential for current conservation, decoupling of spurious states, and for a 



quantitative description of excitation energies of giant resonances [39| 

The residual interaction is derived from the DD-ME2 effective Lagrangian 4o| . 
and single-particle occupation factors at finite temperature are included in a consis- 
tent way both in the FTRMF and FTRRPA. The same interaction is used both in 
the FTRMF equations that determine the single-nucleon basis, and in the matrix 
equations of the FTRRPA. The full set of FTRRPA equations is solved by diagonal- 
ization. The result are excitation energies E v , and the corresponding forward- and 
backward-going amplitudes X u and Y u , respectively, that are used to evaluate the 
transition strength for a multipole operator Qj\ 



B(J,E V ) 



(7) 



Discrete spectra are averaged with a Lorentzian distribution of arbitrary width (r = 1 
MeV in the present calculation). 

In Fig. [1] we display the FTRRPA strength distributions in 56 Ni, 90 Zr, 132 Sn, and 
208 Pb for the isoscalar monopole operator, at temperatures T = 0, 1, and 2 MeV. At 
zero temperature the transition strength distributions are clearly dominated by the 
pronounced collective isoscalar giant monopole resonance (ISGMR). Even by increas- 
ing the temperature to T = 2 MeV, the ISGMR gets only slightly modified in 90 Zr, 
132 Sn, and 208 Pb. In the lighter nucleus 56 Ni, where the ISGMR is less collective and 
therefore more sensitive to modifications of the corresponding configuration space 
at finite temperature, at T = 2 MeV the overall strength distribution is lowered by 
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« 1 MeV. It is interesting to note that for all nuclei under consideration, in the 
temperature range T = 1 — 2 MeV additional transition strength appears at energies 
below 10 MeV. These low-energy transitions occur because of thermal unblocking 
of single-particle orbitals close to the Fermi level. The effect of finite temperature 
is especially pronounced for the neutron-rich nucleus 132 Sn, due to transitions from 
thermally populated weakly-bound neutron orbits. The low-energy strength is dom- 
inated by two rather pronounced peaks. The state at 5.45 MeV corresponds to the 
single-neutron transition from thermally unblocked orbitals: 3p3/2 — > 4p 3 / 2 , whereas 
for the state at 7.02 MeV the dominant transition is between the neutron orbitals 
2f 7/2 and 3f 7 / 2 . 

As already emphasized in the introduction, particularly relevant for possible as- 
trophysical applications is a microscopic description of low-energy dipole response in 
hot nuclei. Fig. [2] displays the FTRRPA (DD-ME2) dipole strength distributions in 
Ni isotopes, at temperatures T = 0,1, and 2 MeV. The transition matrix elements 
are evaluated for the isovector dipole operator. At zero temperature the main peak 
of the isovector giant dipole resonance (IVGDR) is located at ~ 18 MeV. For 68 Ni, in 
particular, additional low-energy PDR structure is predicted below 10 MeV even at 
zero temperature, and this result was recently confirmed in the experimental study 
reported in Ref. |6j. The IVGDR is somewhat modified at finite temperature. At 
T = 2 MeV the main peak is lowered by ~ 2 MeV in 56 Ni, whereas in 62 Ni finite tem- 
perature reduces the fragmentation of the IVGDR strength in the interval ^17-20 
MeV, and enhances the collectivity of the main resonance peak. 

A very interesting finite temperature effect is predicted for the dipole strength 
distributions of 60 ' 62 Ni, where no low-energy excitations are present at zero temper- 
ature. Already at T — 1 MeV novel transition strength develops below 10 MeV. 
For 60 Ni, in particular, at T = 2 MeV the main low-energy peak is located at 
9.71 MeV, it exhausts 1.54% of the Thomas Reiche Kuhn(TRK) sum rule, and is 
composed of 5 neutron and 4 proton single-particle transitions: v2p 3 j 2 — > v2d^/ 2 
(45.87%), ulf s/2 -> v2d 3j2 (13.92%), v2p 3/2 -> u3s 1/2 (8.43%), ulf 7/2 -> ulg 9/2 
(3.58%), p2 Pi/2 -> u2d 3/2 (1.30%), ir2p 3/2 -> n2d 3/2 (7.89%), %2p 3/2 -> 7r2d 5/2 
(6.32%), 7rl/ B /2 -> 7r2d 3/2 (3.62%), and nlf 7/2 -> irlg 9/2 (1.21%). The percent- 
age in brackets corresponds to the contribution of a particular transition to the total 
FTRRPA amplitude J2mi(^mi~^mi)( n i~ n m) for the state at 9.71 MeV. This rather 
rich RPA structure shows that the new low-energy state accumulates a small degree 
of collectivity due to thermal population of both neutron and proton single-particle 
levels around the Fermi surface. 

In the case of 62 Ni, two pronounced low-energy peaks are calculated at temper- 
ature T = 2 MeV: the states at 9.78 MeV and 10.03 MeV, exhausting 1.2% and 
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1.5% of the TRK sum rule, respectively. It is interesting to note that both states 
are dominated by proton transitions. For the state at 9.78 MeV: i\2p 3 j 2 — ► ^2d 5 / 2 
(48.02%), ttI/5/2 -> vr2rf 3/2 (20.99%), ir2p 3/2 -> vr3s 1/2 (3.12%), z/2p 3/2 - ^2rf 5/2 
(17.06%), i/l/ 5/2 -> v2d 3/2 (4.36%), and ulf 7/2 -> i/l0 9/2 (1.22%). The structure of 
the state at 10.03 MeV appears to be less distributed, i.e. it is dominated by the 
transitions v2p 3 j 2 — > v2d 5 / 2 (26.55%), v2p 3 j 2 — >■ ^3si/ 2 (3.59%), 71I/5/2 — > n2d 3 / 2 
(63.60%), and vr2j9 3/2 -> ^2d 5 /2 (2.55%). 

In 68 Ni, where the PDR structure is present already at T — 0, the increase in 
temperature leads to fragmentation of the PDR and its spreading to even lower 
energy. These examples nicely illustrate the variety of effects one can expect for 
the dipole response of hot nuclei, e.g. the appearance of novel modes of low-energy 
excitations dominated by proton transitions, and the thermal spreading of PDR 
structures in neutron- rich nuclei. 

Finally, Fig. [3] shows the evolution of FTRRPA dipole transition strength with 
temperature for 132 Sn. This is an important example, because for this nucleus the 
occurrence of PDR was predicted in several theoretical studies [l|, and recently con- 
firmed in Coulomb dissociation experiments With the increase of temperature up 
to T ~ 1 MeV the transition strength is virtually unchanged. Higher temperatures, 
however, have a pronounced effect on low-energy dipole strength. At temperatures 
pa 2 MeV, even though the main PDR peak at E pa 8 MeV does not change much, ad- 
ditional low-energy single-particle transitions appear in the interval E pa 2 — 1 MeV: 
E=3.75 MeV, z/3p 3/2 -> u3d 5/2 (99.65%); E=4.54 MeV, u2f 7/2 -> u3d 5/2 (97.49%); 
E=4.92 MeV, ^2/ 5/2 -> u2g 7/2 (99.58%); E=6.12 MeV, u2f 7/2 -> u2g 9/2 (98.35%). 
At zero temperature the PDR state is calculated at E — 7.75 MeV, whereas at 
T = 2 MeV the PDR excitation energy is E = 7.77 MeV. In Table 1 we compare the 
structure of FTRRPA states for the main PDR peaks at temperatures T = MeV 
and T = 2 MeV. The amplitudes of transitions dominant at T = MeV get slightly 
reduced at higher temperature. Because of thermal unblocking of the vlg$/ 2 orbit, at 
T = 2 MeV an additional non-negligible contribution to the PDR peak corresponds 
to the transition ^lgg/2 — > vl\in/ 2 . 

The results for the evolution of dipole response with temperature (c.f. Figs. [2] 
and [3]) show that the principal effect of finite temperature on low-energy strength 
that is already present at zero temperature, e.g. in 68 Ni and 132 Sn, is the spreading 
of this structure to even lower energy and the appearance of states that correspond 
to thermally unblocked transitions. This could be important for calculation of astro- 
physical neutron capture rates relevant for r-process nucleosynthesis, because these 
reactions take place at finite temperature. As recently pointed out in Ref. 10|, neu- 



tron capture rates computed with the Hauser-Feshbach model are sensitive to the 
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T = MeV T = 2 MeV 

E = 7.75 MeV E = 7.77 MeV 

v3s 1 /2 — > 3p 3 /2 51.85 50.37 

z/2rf 3/2 -> 3p 3 /2 19.15 17.22 

z/2d 3/2 -> 3pi /2 11.56 10.43 

i/3si/ 2 -> 3pi /2 6.64 5.48 

vlhxi/2 -> lzi3/ 2 4.99 4.44 

7rl# 9/2 -> l/i 11/2 2.17 2.43 



Tabic 1: FTRRPA transition amplitudes for the main PDR peaks in 132 Sn at temperatures T = 
MeV and T — 2 MeV. Included are the contributions (in %) of dominant configuration to the total 
sum of FTRRPA amplitudes: Emi( X mi ~ Y L)i n i ~ 



fine structure of low-energy dipole transitions, emphasizing the importance of reli- 
able microscopic predictions for the evolution of PDR structures both with neutron 
number and nuclear temperature. 

This work was supported by the Unity through Knowledge Fund (UKF Grant 
No. 17/08) and MZOS - project 1191005-1010. Y. F. Niu acknowledges support from 
the Croatian National Foundation for Science, Higher Education and Technological 
Development. The work of J.M and D.V. was supported in part by the Chinese- 
Croatian project "Nuclear structure far from stability" . 

References 

[1] N. Paar, D. Vretenar, E. Khan, and G. Colo, Rep. Prog. Phys. 70 (2007) 691 

[2] D. Savran, M. Babilon, A. M. van den Berg, M. N. Harakeh, J. Hasper, A. 
Matic, H. J. Wortche, and A. Zilges, Phys. Rev. Lett. 97 (2006) 172502 

[3] R. Schwengner, G. Rusev, N. Tsoneva, et al, Phys. Rev. C 78 (2008) 064314 

[4] B. Ozel, J. Enders, H. Lenske, et aI.. larXiv:0901.2443l rl 

[5] P. Adrich, A. Klimkiewicz, M. Fallot, et al., Phys. Rev. Lett. 95 (2005) 132501 

[6] O. Wieland, A. Bracco, F. Camera, et. al., Phys. Rev. Lett. 102 (2009) 092502 

[7] J. Piekarewicz, Phys. Rev. C 73 (2006) 044325 



7 



[8] A. Klimkiewicz et al., Phys. Rev. C 76 (2007) 051603(R) 

[9] S. Goriely, Phys. Lett. B 436 (1998) 10 
[10] E. Litvinova et al., Nucl. Phys. A 823 (2009) 26 
[11] S. Goriely E. Khan, Nucl. Phys. A 706 (2002) 217 

[12] O. Civitarese, R. A. Broglia, and C. H. Dasso, Ann. Phys. 156 (1984) 142 

[13] D. Vautherin and N. Vinh Mau, Nucl. Phys. A 422 (1984) 140 

[14] W. Besold, P.-G. Reinhard, and C. Toepffer, Nucl. Phys. A 431 (1984) 1 

[15] M. E. Faber, J. L. Egido, and P. Ring, Phys. Lett. B 127 (1983) 5 

[16] P. Ring, L. M. Robledo, J. L. Egido, and M. Faber, Nucl. Phys. A 419 (1984) 
261 

[17] H. Sagawa and G. F. Bertsch, Phys. Lett. B 146 (1984) 138 

[18] D. Lacroix, P. Chomaz, and S. Ayik, Phys. Rev. C 58 (1998) 2154 

[19] M. Gallardo, M. Diebel, T. Dossing, and R. A. Broglia, Nucl. Phys. A 443 (1985) 
415 

[20] P. F. Bortignon, R. A. Broglia, G. F. Bertsch, and J. Pacheco, Nucl. Phys. A 
460 (1986) 149 

[21] N. D. Dang, Phys. Rep. 264 (1996) 123 

[22] W. E. Ormand, P. F. Bortignon, and R. A. Broglia, Phys. Rev. Lett. 77 (1996) 
607 

[23] Nguyen Dinh Dang and Fumihiko Sakata, Phys. Rev. C 55 (1997) 2872 

[24] A. B. Larionov, M. Cabibbo, V Baran, and M. Di Toro, Nucl. Phys. A 648 
(1999) 157 

[25] A. N. Storozhenko, A. I. Vdovin, A. Ventura, and A. I. Blokhin, Phys. Rev. C 
69 (2004) 064320 

[26] J. Meyer, P. Quentin, and M. Brack, PHys. Lett. B 133 (1983) 279 
[27] M. Barranco, A. Polls, and J. Martorell, Nucl. Phys. A 444 (1985) 445 

8 



[28] A. L. Goodman, Nucl. Phys. A 352 (1981) 30 

[29] H. M. Sommermann, Ann. Phys. 151 (1983) 163 

[30] E. Khan, Nguyen Van Giai, and M. Grasso, Nucl. Phys. A 731 (2004) 311 

[31] A. Bracco et al., Phys. Rev. Lett. 62 (1989) 2080 

[32] T. S. Tveter et al., Phys. Rev. Lett. 76 (1996) 1035 

[33] D. Pierroutsakou et al, Nucl. Phys. A 600 (1996) 131 

[34] T. Baumann et al., Nucl. Phys. A 635 (1998) 428 

[35] Y. K. Gambhir, J. P. Maharana, G. A. Lalazissis, C. P. Panos, and P. Ring, 
Phys. Rev. C 62 (2000) 054610 

[36] G.A. Lalazissis, J. Konig, and P. Ring, Phys. Rev. C 55, (1997) 540 

[37] T. Niksic, D. Vretenar, P. Finelli, P. Ring, Phys. Rev. C 66 (2002) 024306 

[38] P. Bonche, S. Levit, and D. Vautherin, Nucl. Phys. A 427 (1984) 278 

[39] N. Paar, P. Ring, T. Niksic, and D. Vretenar, Phys. Rev. C 67 (2003) 034312 

[40] G. A. Lalazissis, T. Niksic, D. Vretenar, and P. Ring, Phys. Rev. C 71 (2005) 
024312 



9 







10 



20 

E (MeV) 



30 



40 



10 



20 

E (MeV) 



30 



40 



Figure 1: Isoscalar monopole transition strength distributions in 56 Ni, 90 Zr, 132 Sn, and 208 Pb, 
calculated with the FTRRPA (DD-ME2) at temperatures T = 0, 1, and 2 MeV. 
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Figure 2: Isovector dipole strength distributions in Ni isotopes at temperatures T — 0,1, and 2 
MeV, calculated with the FTRRPA (DD-ME2). 
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Figure 3: Evolution of electric dipole transition strength with temperature in 132 Sn, calculated with 
the FTRRPA (DD-ME2). 
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